Observation of the Nuclear Magnetic Octupole Moment of 87 Rb from Spectroscopic 

Measurements of Hyperfine Intervals 
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The magnetic octupole moment of 87 Rb is determined from hyperfine intervals in the 5p 2 P 3 /2 
state measured by Ye et al. [Opt. Lett. 21, 1280 (1996)]. Hyperfine constants A = 84.7189(22) MHz, 
B = 12.4942(43) MHz, and C = -0.12(09) kHz are obtained from the published measurements. 
The existence of a significant value for C indicates the presence of a nuclear magnetic octupole 
moment Q. Combining the hyperfine constants with atomic structure calculations, we obtain Q = 
—0.58(39) bfiN- Second-order corrections arising from interaction with the nearby 5p 2 P\/2 state 
are found to be insignificant. 
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I. INTRODUCTION 



We report values of the nuclear magnetic octupole mo- 
ment of 87 Rb determined from measured hyperfine in- 
tervals. During an earlierjnvestigation of the hyperfine 
structure of atomic 133 Cs [l[, we discovered existing mea- 
surements of the hyperfine structure of the 5p 2 P 3 /2 state 
of 87 Rb (I = 3/2, F = 0,1, 2, 3) by Ye et al. 0] with 

The measured 
i are shown in 



a frequency resolution of several kHz. 
hyperfine intervals AWf = Wf — Wf 
Fig. [T] These intervals were decomposed in terms of nu- 
clear magnetic dipole and electric quadrupole coupling 
coefficients (A and B) in Ref. Q and the resulting val- 
ues for A and B are given in the first row of Table Q] 
However, for an atomic state with angular momentum 
J > 3/2, there also exists a coupling between the nuclear 
magnetic octupole moment and the electronic state, pro- 
vided I > 3/2. Precision measurements of the three hy- 
perfine intervals in 87 Rb are sufficient to determine three 
hyperfine constants A, B and C, representing the inter- 
action of the nuclear magnetic dipole moment /i/ , electric 
quadrupole moment Q, and magnetic octupole moment 



5p 2 P 3/ 2 




F = 3 



266.650(9) MHz 



F = 2 



156.947(7) MHz 
72.218(4) MHz" 



F = 1 



F = 

FIG. 1: Hyperfine structure of the 5p 2 Ps/2 state in 87 Rb from 
Ye et al. 121 . 
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f2, respectively, with the atomic electrons. We combine 
the coefficient C, determined herein, with atomic struc- 
ture calculations for neutral 87 Rb to determine the nu- 
clear octupole moment. The purpose of the present paper 
is to determine the previously unknown nuclear magnetic 
octupole m oment of 87 Rb, not to discuss the experimen- 
tal work by lYe et all 

The (unperturbed) hyperfine structure intervals of an 
atom having a single electron outside closed shells and 
/ = J = 3/2 is given in terms of the hyperfine constants 
A, B and C by Q] 



AW 3 = 3A + B + 8C 
AW 2 = 2A-B- 28C 
AW! = A-B + 56C. 



(1) 
(2) 
(3) 



Inverting these equations leads to the values of A, B and 
C given in the second row of Table |TJ 



II. SEMI-EMPIRICAL ANALYSIS 

With the aid of approximations developed by Casimir 
semi-empirical estimates for the hyperfine coupling 
constants can be expressed in terms of the nuclear mo- 
ments hi, Q, and fl of 87 Rb. To evaluate the hyperfine 
constants empirically, we start from the nonrelativistic 
expressions for A, B and C @, Chap. VIII], which can be 
written 



A(5P3/2 ) = 



16 / 1 



45 \r 



F x 95.4107 MHz (4) 



Op 



B(5p 3/2 ) = Q(b) /4j ) H 2>> -I. '.)()■) MHz (.V) 



C(5p 3/2 ) = ^(mjv&) 



35 



5p 

R%{r) 



T x 3.40718 Hz. (6) 



The atomic matrix elements in Eqs. (|4|5p are expressed 
in atomic units. In Eq.©, Rbp{ r ) is the radial wave 
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TABLE I: Hyperfine coupling constants A, B and C from 
measured hyperfine intervals of 87 Rb 5p 2 P%/2- First row: re- 
sults obtained by Ye et al. 0, second row: present analysis, 
third row: semi-empirical estimates, fourth row: MBPT re- 
sults. In the semi-empirical and MBPT calculations, we use 
Hi = 2.7516pjv from Raghavan [f| and Q — 133.5(5) mb from 
Pyykko 7]. Q is expressed in ^tjvb. 





A (MHz) 


B (MHz) 


C (kHz) 


Ref. [2] 


84.7185(20) 


12.4965(37) 




Present Work 


84.7189(22) 


12.4942(43) 


-0.12(09) 


Semi-empirical 


91.4 


12.5 


0.293 Q 


MBPT 


83.0 


12.6 


0.206 n 



function of the bp state and, again, atomic quantities are 
expressed in atomic units. The factors F, R, T in the 
above equation s account for relativistic effects. With aid 
of the ICasimirl approximations, we write 









w 


lip 


3(n*) 3 


KM 




4ZfZ 3 


r 4 





9(n*) 3 



_^ da 

dn 

^ da 

dn 



= 0.95983 flg 3 (7) 
= 1393. 68 a^ 5 . (8) 



In the above expressions, Z = 37 — 4 = 33 is the effective 
nuclear charge [3|, p. 166], Z{ = 1 is the ionic charge, n* 
is the effective quantum number of the 5p 3 / 2 state, and 
a is the quantum defect. From Van Wijngaarden ||, 
we find n* = 2.293 and da/dn = —0.052, for the 5p 3 / 2 
state of Rb. Using the values /ij = 2.7516^at [|| and 
Q = 0.1335 b for 87 Rb, we obtain 



A(5p 3/2 ) ~ 89.6 F MHz 
B(5p 3/2 ) » 12.0 i? MHz 
C(5p 3/2 ) « 0.271 T n kHz. 



(9) 
(10) 
(11) 



We estimate the relativistic correction factors F, R, T by 
comparing relativistic and nonrelativistic Hartree-Fock 
(HF) calculations of A, B and C. In this way, we obtain 
F = 1.02, R = 1.04, and T = 1.08. The semi- empirical 
values listed in Table Q] include these relativistic correc- 
tions. As seen from the table, the semi-empirical value 
for A differs from experiment by 8%, while the value of 
B differs by less than 0.1%. 



III. RELATIVISTIC MANY-BODY ANALYSIS 

In Ref. [8[, a relativistic many-body method, referred 
to as the SD approximation, in which single and double 
excitations of Dirac-Hartree-Fock wave functions are it- 
erated to all orders in perturbation theory, was used to 
predict magnetic-dipole hyperfine constants for low-lying 
states of alkali- metal atoms to within a few percent. The 
SD value for the magnetic-dipole hyperfine constant of 
the 5p 3/2 state of 85 Rb from || is A = 24.5 MHz. This 



value is scaled by the ratio 

ff/ [ 87 Rb] _ 2.7516/(3/2) 
ff/ [ 85 Rb] ~ 1.3534/(5/2) 



3.3885 



giving A — 83.0 MHz for the 5p 3 / 2 state of 87 Rb shown in 
the fourth row of TableQ] The all-order method described 
in H was used here to obtain B/Q = 94.16 MHz and 
C/Cl = 0.206 kHz, leading to the values of B and C listed 
in the fourth row of Table |TJ The theoretical uncertainty 
in these values is estimated to be 2%. 

Experimental values of the nuclear moments are listed 
in the first row of Table [TXJ By comparing the all-order 
SD values of Aj [ii and B/Q with the experimental values 
of A, B and C shown in the second-row of Table QJ we 
obtain the MBPT values of nuclear moments fj,i, Q given 
in the second row of Table [TTJ The values of /xj and Q 
determined in this way agree to within 2% with precise 
measurements. Comparing the theoretical value of C/fi 
with the experimental value of C given in the second row 
of Table Q] leads to the principal prediction of the present 
paper: Q = -0.58(39) ^wb. 



IV. ESTIMATING NUCLEAR MOMENTS 

It is of interest to compare values of nuclear moments 
inferred from atomic structure calculations with values 
obtained directly from nuclear shell- model calculations. 
In the extreme shell model [9J, properties of 87 Rb can 
be described assuming a single valence nucleon moving 
around an inert core. According to Schwartz (loj . the 
shell- model predictions for the nuclear moments fii, Q 
and Q are: 



M = fj, N Ix 

' 9L + {gs~9L)/{2I), 



I = L + l/2 



9L -(g s -giL)/(2I+l), I = L-l/2 



Q = -^—^9L{r 2 ) 
^ 2I + 2 y X 1 



_ 2 (27-1). 2 . 
3' (2/ + 4) V 1 

(I + 2)[(I-3/2)g L + g s ], I = L + l/2 
{I-mi + h/2)g L - 9s ], I = L-\/2 



(12) 



(13) 



(14) 



For 87 Rb, the unpaired p 3 / 2 proton has total angular mo- 
mentum I = 3/2, orbital angular momentum L = 1, and 
spin angular momentum S — 1/2. The proton spin gy- 
romagnetic ratio is gs = 5.585694701(56) [ll[ and the 
orbital gyromagnetic ratio is g^ = 1. Using the value 
of the mean-squared nuclear radius (r 2 ) = 0.180 b from 
[l2| , we obtain the values of the three nuclear moments 
given in the third row of Table [Til Comparing the shell- 
model value of Q with the experimental value shown in 
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TABLE II: Nuclear moments of 87 Rb. First row: mean exper- 
imental values, second row: values calculated using MBPT 
matrix elements and the hyperfine constants given in Ta- 
ble HI third row: values calculated in the extreme nuclear 
shell model Eqs. (|12ll4|) . 





in (maO 


Q (mb) 


Q (fj, N b) 


Expt. 


2.751639(2) 


133.5(5) 




MBPT 


2.70 


135 


-0.58(39) 


Shell-model 


3.79 


-72 


0.30 



Table HT1 it is clear that nuclear many-body effects modify 
both sign and magnitude of the shell-model prediction. 
Therefore, the difference in sign and magnitude between 
the value of £1 predicted here and given in the second row 
of Table Hi] and the shell-model value shown in the third 
row is not particularly surprising! 



V. SECOND-ORDER CORRECTIONS 

Second-order hyperfine corrections will modify Eqs. {1} 
[3]) and possibly influence the value of C extracted from 
these equations. In this section, we investigate the in- 
fluence of second-order corrections and show that they 
have a negligible effect on the 5p 3 / 2 hyperfine constants 
in 87 Rb. The second-order correction to the energy of a 



state |0) can be written 



|H hf |n)(n|ff hf |0) 
Eq — E„ 



(15) 



where H^f is the hyperfine interaction Hamiltonian. For 
the case of interest here, |0) is a particular hyperfine sub- 
state of the 5^3/2 state and the sum over states |n) is re- 
stricted to substates of the nearby 5pi/ 2 state. Only the 
two substates F = 1 and F = 2 of the 5p 3 / 2 states are 
modified by interaction with the 5pi/ 2 state. Following 
the discussion in [131, we find 



W, 



(2) 



(16) 



where 



2 20 <5p 3 /2||Ti||5p 1/2 ) 

v = "'T ~ae 



(17) 



20V3 (5p3/2\\Ti\\5 Pl/2 )(5p 3/2 \\T 2 \\5p 1/2 



AE 



(18) 



The magnetic dipole operator T\ and electric quadrupole 
operator T 2 in the atomic reduced matrix elements in 
the numerators of Eqs. (fTT)) and (fT5|) are defined in [14j, 
Chap. V], and the energy denominator AE in Eqs. (fTT|) 



and (|18p is the 5pz/ 2 — 5pi/ 2 fine-structure interval. Rel- 
ativistic many-body calculations carried out in the SD 
approximation [H give 

(5p 3 /2||7i||5p 1/2 ) = 22.3MHZ//MT 
(5p 3 /2l|T 1 2||5p 1/2 ) = 219.6 MHz/b. 

The 5p fine-structure interval in 87 Rb is AE = 
7,123,020.80(5) MHz [Hj]. Combining these values, we 
find i] = 3.524 kHz and C = 2.916 kHz. 

Including second-order corrections, Eqs. become 



AW 3 = 3A + B 



20 ' 100 



AW 2 



2A- B 



28C+1,- 



AWi = A - B + 56C 



36 



2^/5 
" ~75~ ( 
V5 
60 



c 



(19) 
(20) 
(21) 



Inverting these equations leads to the following second- 
order corrections to the previously determined values of 
A, B and C listed in the second row of Table U 



A (2) 
C {2) 



180 



n - 



V5 

750 



C = 0.0000109 MHz 



—n + — C = 0.000183 MHz 
30 ' 100 s 



V5 
2000 



C = 0.00326 kHz. 



(22) 
(23) 
(24) 



The above equations for the second-order corrections are 
in precise agreement w ith results for 87 Rb presented by 
iBelov and Dereviankol in Ref . [l6| . The second-order cor- 
rections leave the experimental values of A and C listed in 
the second row of Table |T] unchanged and insignificantly 
increase the value of B to 12.4944(43) MHz. 



VI. CONCLUSIONS 

In summary, the precise hyperfine measurements of the 
5p 3 / 2 hyperfine interval in 87 Rb by Ye et al. are rean- 
alyzed, allowing for the possibility of a nonvanishing nu- 
clear magnetic octupole moment. Assuming that second- 
order effects in the hyperfine interaction are negligible, 
we obtain the values listed in the second row for the hy- 
perfine constants A, B and C. In particular, we find a 
nonzero value for C = —0.12(09) MHz. Values of A/ni, 
B/Q and C/Q are evaluated both empirically and using 
relativistic all-order methods. The empirical calculations 
in combination with the experimental measurements lead 
to values of \ii and Q that agree with other measured 
values to better than 8%, while the all-order calculations 
combined with the experimental hyperfine constants lead 
to values of /ij and Q that agree with other measurements 
at the 2% level. With the aid of all-order calculations, 
we infer the value fl = —0.58(39) /i„b for the nuclear 
magnetic octupole moment. This value (together with 



4 



the measured value of Q) is larger in magnitude and dif- 
ferent in sign than the value predicted by the nuclear 
shell model. We examined the influence of the second- 
order hypcrfine corrections arising from the interaction 
between the 5p 3 / 2 state and the neighboring 5pi/ 2 state. 
These corrections are found to make insignificant changes 



in the values of A, B and C extracted from the measure- 
ments. Inasmuch as the experimental uncertainties in 
the values of C and f2 are relatively large, more precise 
measurements of the hyperfine intervals would certainly 
be desirable. 
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